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Abstract 

We study the entanglement of formation for arbitrary dimensional bi- 
partite mixed unknown states. Experimentally measurable lower and 
upper bounds for entanglement of formation are derived. 

PACS numbers: 03.67.-a, 02.20.Hj, 03.65.-w 

Being one of the most striking phenomena in quantum physics, quantum entanglement 
[l-3] has been extensively investigated in recent years. One of the main tasks in quantum 
entanglement theory is to quantify the entanglement of quantum systems. Among all the 
bipartite entanglement measures, entanglement of formation (EOF) is one of the most mean- 
ingful and physically motivated measures 0, [sj , which quantifies the minimal cost needed to 
prepare a certain quantum state in terms of EPR pairs, and plays important roles in many 
physical systems, such as quantum phase transition for various interacting quantum many- 
body systems {fjj], macroscopic properties of solids (tJ, and capacity of quantum channels 

3. 

Let Ha, Hb be the m, n (m < n) dimensional vector spaces respectively. A pure quantum 
state G Ha ® Hb is an mn- dimensional vector. Its entanglement of formation is defined 
by E(\ip)) = S(pa), where pa = TrB(\ip)(ip\) is the reduced density matrix of S(pa) 
is the entropy 

m 

S(pa) = -^Pihgpi = H(p), (1) 

i=i 

where log stands for the natural logarithm throughout the paper, p^ are the eigenvalues of pA 
and p is the Schmidt vector (pi, p 2 , ■ ■ ■ ,p m )- This definition of entanglement of formation 
is extended to mixed states p by the convex roof, 

E(j>)= mm ^ptEdin)), (2) 
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for all possible ensemble realizations p = J2iPi\i J i)(4'i\i where Pi > and YliiPi = 1- 

It is a great challenge to find an analytical formula of the entanglement of formation 
for general bipartite quantum mixed states p. Considerable efforts have been made on 
deriving entanglement of formation or its lower bound through analytical and numerical 
approaches. So far the entanglement of formation has been calculated for some particular 
states like bipartite qubit states [9j, isotropic states [l0| and Werner states [ll| in arbitrary 



dimensions, and symmetric gaussian states in infinite dimensions 12j. In 13J in order to 



estimate the entanglement of formation for general states, a lower bound has been presented 
by using the partial transposition of p with respect to the subsystem H A , p TA , and the 
realignment of p, R(p). It is shown that (without regard to the normalization coefficient 
log 2) 

0, n = 1; 

ff 2 [7(«)] + [1 - 7(ft)] log 2 (™ - 1), K < 



E(P) > < 



4(m— 1) . 



(3) 



log 2 (m-l) 
m-2 



(fi — m) + log 2 m, 



4(m-l) 



< < m; 



where $7 = max{||p TA ||, ||P(p)||}, H 2 is the standard binary entropy function, \\A\\ denotes 
the trace norm of the matrix A. 

Comparing with the entanglement of formation, the entanglement measure concurrence 
is relatively easier to be dealt with. In [14] a simpler analytical lower bound for concurrence 



has been presented. And a series of new resu 



been further obtained [15( . In particular in 
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ts related to the bounds of concurrence have 



3Ml8l|. the authors derive measurable lower and 
upper bounds for concurrence for general mixed quantum states, 

2[Trp 2 - Trp^] = Tr(p g> P V$ < [C{p)f < Tr(p <g> pK { ) = 2[1 - Trp 2 A ], (4) 

where V x = 4(P_ -P+) ®P_, V 2 = 4P_ ® (P_ - P+) and K x = 4P_ ® I, K 2 = 4/(8) P_, P_ is 
the projector on the antisymmetric subspace of the two copies of either subsystem, P + the 
symmetric counterpart of P_. 

Contrary to the concurrence, less has been achieved related to the lower and upper bounds 
of entanglement of formation. In this paper, we derive analytical lower and upper bounds 
for entanglement of formation which care measurable experimentally. These bounds supply 
nice estimation for entanglement of formation for some quantum states. 

To derive lower and upper bounds for entanglement of formation, we first consider a pure 
|^) with Schmidt decomposition \ip) = \/JhV^)^ where p« > 0, Y1T fa = 1- ^ ^ s easn y 
verified that 

1 - Trp^ = Trp 2 - Trp^ = 1 - J> 2 = A. (5) 



Set 



X(X) = max{P(/i)|l = ^(A) = mm{H(p)\l -J2^ = A >- 



(6) 
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Let e(x) be the largest convex function that is bounded above by Y(x) and r](x) the smallest 
concave function that is bounded below by X(x). 

Theorem: For any m®n(m < n) quantum state p, the entanglement of formation E(p) 
satisfies 

max{e(Trp 2 - Trp 2 A ), e{Trp 2 - Trp|)} < E(p) < min{r/(l - Trp 2 A ), 77(1 - Trp|)}. (7) 

Proof: Without lose of generality, we assume that e(Tip 2 — Trp A ) > e(Trp 2 — Trp 2 ^) and 
77(1 — Trp 2 ^) < f](l — Trp^). Note that for any pure state the concurrence is given by 



C (1^) ) = y2(l — Tr(p A ) 2 ). Due to convexity of concurrence, for any pure decomposition 
P = J2 a PaPa, we have 

a a 

Taking into account of the bounds of C 2 (p), we obtain 

5>«[Tty>«) 2 - Tr(p^) 2 ] > Trp 2 - Trp 2 A (8) 

a 

and 

E^t 1 -^) 2 ]^!-^- (9) 

a 

Assume p = J2 a PaPa be the optimal decomposition of E(p). We have 

a a a a 

where A = Trp 2 — Trp 2 ^. We have used the definition of e to obtain the first inequality. The 
second inequality is due to the convex property of e(x) and the last one is derived from (jSJ). 
On the other hand, 

E{ P ) = E^dvo) = ^ Y.p^(K) < »k$>«a«) ^ ^ A ')> ( n ) 

a a a a 

where A' = l — Tip A . We have used the definition of r\ to get the first inequality. The second 
inequality is derived from the concave property of rj(x) and the last one is obtained from 
©. □ 
We calculate now both the maximal admissible H(p) and the minimal admissible H(p) 



for a given A, i.e. X(X) and Y(X), by using the Lagrange multipliers approach lOj. The 
necessary conditions for the maximum and minimum are given by: 

- log p k - 1 - 2xp k + y = 0, (12) 
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1-J2tf-A = 0, 1-^ = 0; (13) 

i i 

where x, y denote the Lagrange multipliers. We get from (fT2|) that 

log^ fc = -2xjj, k +y - 1. (14) 

From (I14p we know that there are at most two solutions for each which will be denoted 
as a and 8 in the following. 

Let rii be the number of entries where fii = a and n 2 the number of entries where /ij = 8. 

If one of the ni and no is zero, we have that A = 1 7 — . Otherwise both a and 8 are 

nonzero. The problem is now turned to be, for fixed n 1; n 2 , n\ + n 2 < m, one maximizes or 
minimizes the function 

(A) = mhipt) + n 2 h(8), (15) 
where h(x) = —xlogx, under the constraints (ITS]) . By direct computation we obtain 



± 



± _ wi ± y^j - 7ii(wi + w 2 )[l - n 2 (l - A)] ± _ 1 - nig , . t 

ni(ni + n 2 ) n 2 

To ensure the nonnegativity property of a^ in2 and 8^ lTl2 , we require that max{l , 1 — 

±} < A < 1 - j^i^. Since a~ ni = ^ ina ,/3~ ni = a+ n2 , the function in Eq. $T5]) takes 
the same value for ai n2 and a;~ x . Therefore we can restrict ourselves to the solutions 
d ni n 2 — a n in2 - ^q. ([15T) then turns out to be 

K in2 W = nih(a+ im ) + n 2 h(8+ in2 ). (17) 

When m — 3, to find the expressions of upper and lower bounds in ([7]) is to obtain the 
max- and minimization over the three functions iq 2 (A), F 2 i(A) and Fn(A). From (ITTj) for 
m = 3 we have 

(F 12 , i<A<| [F 21 , 1<A<|. 

From (ITS]) we have that 77 [A] is the broken line connecting the following points: 
[0,0],[i,log2],[|,log3]. 

In order to determine e[A] we solve the following equations: Let 1(A) = k(A — 0.5) +0.868 
be the line crossing through the point [0.5, iq 2 (0.5)]. We solve (i) /(A) = Fu and (ii) 
^dlT~ = k = dF ^ A ' > for k and A, and find the values to be 1.65 and 0.091. Thus we derive 
that e[A] is the curve consisted of Fu for < A < 0.091 and a broken line connecting points 
[0.091, Fn(0.091)], [0.5, F 12 (0.5)] and [0.667, log[3]], i.e. 

, 21og2 x A, < A < 0.5; 
77 [A] = < " (19) 

' 6 log § x (A — |) + log 2, 0.5<A<0.667 
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FIG. 1: Upper and lower bounds of E{p) (solid lines) for m = 3, and F\\, F12, F21 (dashed line). 



and 

F u , < A < 0.091; 
e[A] = 1 1.65(A - 0.5) + 0.868, 0.091 < A < 0.5; (20) 
1.39(A - 0.667) + 1.099, 0.5 < A < 0.667, 

see Fig{TJ 

When m = 4, we need to find the max- and minimization over the six functions 
Fu, F12, F21, F22, F31, and F13, which are plotted in Figj2j We have 



X(A) 



F n , 0<A<i 

F\2 
F13 



2 ' 

|<A<1; 
1<A<1 



y(A) 



F 1U 0<A<| 
F21 
F 



2 ■ 

2 ^ 1V — 3' 



(21) 



31, I < A< f. 



Further more, one obtains that 77 [A] is the broken line connecting the following points: 
[0,0], [|,log2], [|,log3], [§,log4] and e[A] is the curve consisted of F u for < A < 0.062 and 
a broken line connecting points [0.062, 0.142], [0.667, 1.242] and [0.75, 1.386] (These points can 
be derived by using the same processes as that have been done in case m — 3), i.e. 



V [A] 



21og2xA, < A < 0.5; 

6 log § x (A- \) + log 2, 0.5<A< 0.667; 



(22) 



12 log I x (A 



log 3, 0.667 <A<0.75 



and 



e[A] 



Fu, < A < 0.062; 

1.820(A - 0.667) + 1.242, 0.062 < A < 0.667; 
1.726(A - 0.667) + 1.242, 0.667 < A < 0.75. 



(23) 
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FIG. 2: Upper and lower bounds of E(p) (solid lines) for m = 4, and i*ii, -F12, -^21, 
^22, ^13,-^31 (dashed line). 

Generally, we have the following observation, for any m, 



^11, 0<A<i; 

*i2, i<a<|; 



X(A) = { 



Y(A) = 



F u , 0<A<±; 
F21, I < A< |; 



k Fi {m -i), < A < 



m— 1 . 

m ! 



(24) 



77 m— 2 ^ \ ^ m—1. 

I ^(m-l)l, T^TT < A < — , 



and 77 [A] is the broken line connecting the following points: [-^-, log (i + 1)], < % < m — 1, 
i.e. 

r/[A] = fc(fc + 1) log *±I(A - ^) + log fc, (25) 

for (k — 1) < A < k and = 1, 2, • • • , m — 1. The representation of s[A] can be also figured 
out numerically. 

The measurable upper and lower bounds can be used to estimate the entanglement of 
formation for an unknown quantum mixed state experimentally. Consider the following 
mixed quantum state 

p=Zl + (l-x)\4>)(4>\. (26) 



where = (a, 0, 0, 0, 0, 0, 0, ^Y/y/Tr{\i/))(i/;\}. For x = 0.1, one has 
\ I 1^1-TrVil-Tff, ! -Tr{,£} = 



a T ,,..2, , .,-,,.,.21 1.45 + 9.21a 2 -0.38a 4 

' " ' ~ (2 + 3a 2 ) 2 



A' = 1 - Tr{pi} = 1 - Tr{p 2 } = 



1.14(0.19 + a 2 )(9.67 + a 2 ) 
(2 + 3a 2 ) 2 ' 



(27) 
(28) 



6 



Substituting A and A' above into (fl9]) and (1201) respectively, we have the upper and lower 
bounds (j7|) for E(p), see Fig. 3. 
For x = 0.001, one has 



A = Tr{p 2 } - Tr{pi} = Tr{p 2 } - Tr{p|} 



1.99+ 11.97a 2 - 0.004a 4 

(2 + 3a 2 ) 2 



2| i rr2) 0.01(0.17 + a 2 )(999.67 + a 2 ) 

(2 + 3a 2 ) 2 



A' = l-Tr{pi} = l-Tr{p 2 } 



(29) 



(30) 



The corresponding bounds of E(p) is shown in Fig. 4. We see that the lower and upper 
bounds are closer. And the value for E(p) can be estimated more precisely. 
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FIG. 3: Upper and lower bounds of E(p) with x = 0.1. 
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FIG. 4: Upper and lower bounds of E(p) with x = 0.001. 
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We have studied the entanglement of formation for mixed quantum states. We have 
derived upper and lower bounds for entanglement of formation that are experimentally 
measurable. These bounds together can be used to estimate the entanglement of forma- 
tion for arbitrary finite dimensional unknown states according to a few measurements on 
a twofold copy p £g> p of the mixed states. These results supplement further the estimation 
for entanglement of formation, like the case of concurrence for which many lower and upper 
bounds have been already obtained. 
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